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Using a sample of 1.31 billion J/i/) events accumulated with the BESIII detector at the BEPCII 
collider, we report the observation of the decay J/ip —>■ ^ which is the first evidence for a 

doubly Okubo-Zweig-Iizuka suppressed electromagnetic J/i!) decay. A clear structure is observed in 
the K^K~ mass spectrum around 1.02 GeV/c^, which can be attributed to interference between 
J/tp —^ (pn^ and J/tp —>■ K^K~ '7r° decays. Due to this interference, two possible solutions are found. 
The corresponding measured values of the branching fraction of J/p) —>■ (p-K^ are [2.94± 0.16(stat.) ± 
0.16(syst.)] X 10“® and [1.24 ± 0.33(stat.) ± 0.30(syst.)] x 10“^. 

PACS numbers: 13.25.Gv, 14.40.Be 


The discovery of the J/ip played an important role 
in understanding the basic constituents of nature 
and opened a new era in particle physics. Its un¬ 
expected narrow decay width provided insight into 
the study of strong interactions. As its mass is be¬ 
low the charmed meson pair threshold, direct decay 
into charmed mesons is forbidden. Therefore the 
J/ip hadronic decay modes are Okubo-Zweig-Iizuka 
(OZI) 0 suppressed, and the final states are com¬ 
posed only of light hadrons. 

A full investigation oi Jjpi decaying to a vector 
meson (V) and a pseudoscalar meson (P) can pro¬ 
vide rich information about SU(3) flavor symmetry 
and its breaking, probe the quark and gluon con¬ 
tent of the pseudoscalar mesons, and determine the 
electromagnetic amplitudes However the pres¬ 

ence of doubly OZI (DOZI) suppressed processes, 
like the observation oi J/pi radiatively decaying into 
ui(p Si, complicates matters as they do not obey 
quark correlation or satisfy nonet symmetry (treat¬ 
ing SU(3) octets and singlet as a nonet and assum¬ 
ing the coupling constants are the same in the inter¬ 
actions 0,0)- Well established phenomenological 
models 0, Isj have indicated that the DOZI ampli¬ 
tude can have a large impact through interference 
with the singly OZI suppressed amplitude. 

Of interest is the decay J/pi ^ (pn^, which occurs 
via the electromametic DOZI process or by non¬ 
ideal ui — (p mixing 1 , 10 . Recently, using a combi¬ 
nation of a factorization scheme for the strong decays 
and a Vector Meson Dominance (VMD) model for 
electromagnetic decays in J/pi —>■ VP, the branch¬ 
ing fraction oi J/pi —>• 6-k^ has been predicted 
to be around 8 x 10“^ Q, while the best upper 
limit to date comes from the BES collaboration, 
B{J/pi —>■ ‘P’^^ < 6.4 X 10“® at the 90% confidence 
level (C.L.) Q. In this paper, we report the first 
observation of J/pi —>■ cpir^ based on a sample of 
(1.311 ± 0.011) X 10® J/pi events [13, [HI accumu¬ 
lated with the BESIII detector. 

The BESIII detector [l3] is a magnetic spectrom¬ 
eter located at the Beijing Electron Positron Col¬ 


lider (BEPCII), which is a double-ring e+e” collider 
with a design peak luminosity of 10®® cm“®s“® at 
the center of mass (c.m.) energy of 3.773 GeV. The 
cylindrical core of the BESIII detector consists of 
a helium-based main drift chamber (MDC), a plas¬ 
tic scintillator time-of-flight system (TOE), and a 
CsI(Tl) electromagnetic calorimeter (EMC). All of 
them are enclosed in a superconducting solenoidal 
magnet providing a 1.0 T (0.9 T in 2012) magnetic 
field. The solenoid is supported by an octagonal 
flux-return yoke with resistive plate counter muon 
identifier modules interleaved with steel. The accep¬ 
tance for charged particles and photons is 93% of 47r 
solid angle. The charged-particle momentum resolu¬ 
tion is 0.5% at 1 GeV/c, and the specific energy loss 
(dE/dx) resolution is 6%. The EMC measures pho¬ 
ton energies with a resolution of 2.5% (5%) at 1 GeV 
in the barrel (endcaps). The time resolution of the 
TOE is 80 ps in the barrel and 110 ps in the endcaps. 
The BESIII offline software system (BOSS) frame¬ 
work is based on Gaudi [T3]. A GEANT4-based [3 
Monte Carlo (MG) simulation is used to determine 
detection efficiencies and estimate backgrounds. 

For the decay J/pi —>■ (pir^ —^ K'^K~'y'y, a can¬ 
didate event is required to have two charged tracks 
with opposite charge and at least two photons. For 
each charged track, the polar angle in the MDC must 
satisfy |cos0| < 0.93, and the point of closest ap¬ 
proach to the e“'"e“ interaction point must be within 
±10 cm in the beam direction and within 1 cm in 
the plane perpendicular to the beam direction. TOE 
and dE/dx information are combined to give parti¬ 
cle identification (PID) probabilities for tt, K and p 
hypotheses. To identify a track as a kaon, the PID 
probability for the kaon hypothesis must be larger 
than that for the pion hypothesis. 

For each photon, the energy deposited in the EMC 
must be at least 25 MeV for | cos0| < 0.8 or 50 MeV 
for 0.86 < I cos 6*1 < 0.92. To select isolated show¬ 
ers, the angle relative to the nearest charged track 
must be larger than 20°. The timing information of 
the EMC is used to suppress electronic noise and 
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FIG. 1. (a) Scatter plot of M[K+K- ) versus M( 77 ). 
The red solid and blue dotted boxes are the tt*’ signal and 
sideband regions. The red dashed box indicates J/i/> —>■ 
07. (b) K+K- mass spectrum, where the dots with 

error bars are events in the 7r° signal region; the hatched 
histogram are events from the 7r° sidebands; and the 
dashed histogram is MC simulation of J/ijj —>■ 07r'^ with 
arbitrary normalization. 


unrelated energy deposits. Furthermore, a four- 
constraint (4C) kinematic fit is applied to the candi¬ 
date events under the hypothesis, requir¬ 

ing the 4-momentum of the final state to be equal 
to that of the colliding beams. If there are more 
than two photon candidates in an event, the combi¬ 
nation with the smallest xlciK'^K~jj) is retained. 
Events with xic < 30 are selected. 

After the above selection, the scatter plot of 
versus M{xx) (Fig- [D (a)) shows two 
clear clusters corresponding to 07 and 0 ^' and two 
bands corresponding to and K'^K~r ], but 

no evident accumulation of events for 07 r*^. To inves¬ 
tigate the M{K^K~) spectrum of events, 

we select events where the 77 invariant mass is in 
the 7 r° mass region 0.115 < M{-fx) < 0.155 GeV/c^. 
The M{K^K~) distribution for these events is 
shown in Fig.[I](b), where a clear structure around 
the 0 mass is seen. 

Studies were performed using both MC events and 
data to investigate whether the structure around 


TABLE I. Background analysis for the decay J/ij) —>■ 
07r’^. 


Type 

Reactions 

Coherent 

J/0 ^ A'+A'-yr® 

0 peaking 

e+e“ 7/sfl0, J/i> 07r®7r®/077 

TT® peaking 

J/0 —>■ 77c(lS') —>■ yA+A^TT® 

Other 

J/0 7A+A-/77r°A+A-/7r®7r®A+A- 


1.02 GeV/c^ could be background related. We ana¬ 
lyze a MC sample of 1.2 x 10® J/ip inclusive decays, 
in which the known decay modes were generated by 
BesEvtGen [illii with measured branching frac¬ 
tions [13 while unknown decays were generated by 
Lund-Charm [T^. The dominant background events 
are found to be from J/^i —?> A'+AT^tt® with the in¬ 
termediate states decaying into AT^tt® and AT+AT”, 
which is coherent for the decay J/ijj —>■ 07 r®. A 
partial wave analysis, not including J/ip —>■ 07 r® 
but considering the interference of all intermediate 
states, yields a smooth distribution with K'^K~ 
mass below 1.2 GeV/c®. The incoherent background 
can be categorized into three classes as follows. ( 1 ) 
The 0 peaking background: e+e" —>■ xisr4> and 
J/0 —>■ 07 r® 7 r®/ 077 . The former background is 
studied using data taken at energies far from any 
charmonium resonance and the latter ones are stud¬ 
ied by exclusive MC samples. The studies show 
that these background events can be compensated 
by the tt® mass sideband events, which are dehned 
as 0.055 < M( 77 ) < 0.095 GeV/c® and 0.175 < 
< 0.215 GeV/c®. (2) The tt® peaking back¬ 
ground: J/ijj ^ 77c( 15') —>■ yAT+AT^TT®. This back¬ 
ground cannot be taken into account by the tt® mass 
sideband events. From MC simulations, the ratio 
between the number of this background events and 
the number of the coherent background events in 
the TT® mass region is 0.5%. As little is known 
about the possible intermediate states, we neglect 
this background and consider the related systematic 
uncertainty. (3) The non-0 and non- 7 r® background 
are dominated by the decays J/tp —>• yAT+A'”, 
77 r®A'+Ar“ and K'^K~ with various interme¬ 
diate states. MC simulations show they can be sub¬ 
tracted by the tt® mass sideband events. All back¬ 
ground types are summarized in Table HI Through 
the studies above, none of these background events 
produce a structure in the K^K~ mass spectrum. 
In addition, the detection efficiency as a function of 
M{K^ K~), obtained from the MC simulation and 
taking into account the angular distributions [l^ , is 
also smooth over the K'^K~ mass region, with no 
structure in the region of the 0 signal. 
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A possible explanation for the structure in the 
M{K'^K~) spectrum is interference between J/ij) —>■ 
07 r° and other processes with the same final state. 
We have verified this using a statistical hypothe¬ 
sis test [lO, [2l|. In the null hypothesis without 
J/ip —>■ 4>t^^ , a second-order polynomial function, 
defined as Fho = P(rn) = cq + Cim -I- C 2 m^, is used 
in the fit to describe the data after subtraction of 
the 7 r°-sideband events. The positive hypothesis is 
characterized by a two-component function in 

which the model is a coherent sum of a relativistic 
Breit-Wigner resonance and the second-order poly¬ 
nomial function, convoluted with a Gaussian func¬ 
tion G{m, (Tm) to take into account the mass resolu¬ 
tion, am- 



Phi = 


where 


\^/P(rnjJ^{^+A^{m)\'^^{rn)^G{m,am), 

( 1 ) 


Arj,{m) 

with 


y/Re'^^p^{m)pKim) B{p^{m)) B{pK{m)) 
m? — vtlqF imT{m) B{p^{mo)) B{pK{mo)) ’ 

( 2 ) 


r(m) 


/ pKjm) y mp B{pK{m)) ^ 
\PK{mo)J m B{pKimoj) 


Here, m is the K^K invariant mass, mo and 
To are the nominal mass and decay width of the 
</> [i3- P<j>{jn) {pKim)) is the momentum of the (j) 
(AT) in the frame of J/ip (</») with the mass of (j) being 
m. $(m) = Ptp(m)pK(rn) is the phase space factor. 
B{p), defined as B{p) = 1 /is the Blatt- 
Weisskopf penetration form factor with the me¬ 
son radius r being 3 GeV“^. R and 5 represents the 
magnitude and relative phase angle respectively for 
the contribution of the (j) resonance. Omitting the 
convolution, Fh^ can be expanded to be 


P(m) -I- |A,^{m)|^$(m) -f 2yJP(m)^{m)^A^{m) , 

(4) 

where the first term is the non-(/) contribution from 
the decay J/ip ^ K^K~tA‘\ the second term is the (p 
resonance from the decay J/ip ^ 07 r°; and the third 
term is their interference. Here, ^A^{m) denotes 
the real part of A^{m). 

MC simulations show that the K~ mass res¬ 
olutions are essentially the same for J/p) decay¬ 
ing to and prj with 'K^/rj —>■ 77 . We obtain 
cTm = (1.00 ± 0.02) MeV/c^ by performing an un¬ 
binned likelihood fit to the M{K^K~) spectrum 
of J/p —)> (pTj with 0.50 < M( 77 ) < 0.60 GeV/c^, 
shown in the red dashed box in Fig. [T] (a). The 


FIG. 2. Fit to M{K^K~) spectrum after sideband sub¬ 
traction for Solution I (a) and Solution II (b). The red 
dotted curve denotes the p resonance; the blue dashed 
curve is the non-p contribution; the green dot-dashed 
curve represents their interference; and the blue solid 
curve is the sum of them. 


same Breit-Wigner formula convoluted with a Gaus¬ 
sian function is used to describe the p signal, while 
a second-order polynomial is used to describe the 
background. 

After subtracting the incoherent background 
events estimated with 7 r° sidebands, a maximum 
likelihood fit is performed to the M{K'^K~) dis¬ 
tribution under the positive hypothesis. Two solu¬ 
tions with two different phase angles between the 
p resonance and the TiOTi-p contributions are found. 
The final fits, including the individual contributions 
of each components, are illustrated in Fig. [5] (a) 
and Fig. [5] (b), while the signal yields and the rel¬ 
ative phase angles are summarized in Table |H1 In 
Fig. [5] (a) and Fig. [5] (b), the blue dashed curve is 
the non-0 contribution (the first term in Eq. IT]); the 
red dotted curve denotes the p resonance (the sec¬ 
ond term in Eq. [d]) ; the green dot-dashed curve rep¬ 
resents their interference (the third term in Eq. 0]) ; 
and the blue solid curve is the sum of them. The 
signal yield AI®*® in Table |TT] is calculated by inte¬ 
grating the function of the p resonance over the fit 
range. The statistical significance is determined by 
the change of the log likelihood value and the num¬ 
ber of degrees of freedom in the fit with and without 
the p signal [13, Both solutions have a statistical 
significance of 6.4tT, which means that they provide 
identically good descriptions of data. 

The non-resonant quantum electrodynamics 
(QED) contribution is estimated in two ways. One 
way is by analyzing data taken at energies far from 
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TABLE II. Fit results. is the fitted number of signal 
events (from the parameter R). <5 is the relative phase. 
2 A log £ is 2 times the difference of the log-likelihood 
value with and without (j) signal, while Nf is the change 
of the number of degrees of freedom. Z is the statistical 
significance. 


Solution 

N-is 

5 

2A log C/Nf 

Z 

I 

838.5 ±45.8 

-95.9° ± 1.5° 

45.8/2 

6.4(7 

II 

35.3 ±9.3 

-152.1° ± 7.7° 

45.8/2 

6.4(7 


any resonance, namely at 3.05, 3.06, 3.08, 3.083, 
3.090 and 3.65 GeV. The other way is to use data 
from the ^/>(3770) resonance, assuming that the 
possible contribution tj}{3770) —>■ 4>Tr^ is negligible. 
The selection criteria are the same except for the 
required 4-momenta in the kinematic fit. Neither 
sample shows significant events. With a 

simultaneous fit, we obtain the QED contribution 
to the signal yield (3.097) < 5.8 at the 90% 
C.L., normalized according to the luminosity and 
efficiency and assuming the cross section is propor¬ 
tional to 1 /s with s being the square of the c.m. 
energy. Thus we neglect the non-resonant QED 
contribution and use the upper limit of (3.097) 
to estimate a systematic uncertainty from this 
assumption. 

With the detection efficiency, (45.1 ± 0.2)%, ob¬ 
tained from the MC simulation, the branching frac¬ 
tions of J/ij} —>■ (/) 7 r° are calculated to be (2.94 ± 
0.16) X 10“® for Solution I and (1.24 ±0.33) x 10“^ 
for Solution II, where the errors are statistical only. 

The sources of systematic uncertainty and their 
corresponding contributions to the measurement of 
the branching fraction are summarized in Table IIIII 
The tracking efhciency of charged kaons is stud¬ 
ied using a high-purity control sample J/'ij} —>■ 
, while the photon detection efficiency is 
investigated based on a clean sample of J/ip —>■ p 7 r. 
The differences between data and MC simulation are 
1 % for each charged track and 1 % for each pho¬ 
ton. The 7 r° selection efficiency is studied with the 
sample J/ij) —>■ pn, and MC simulation agrees with 
data within 0.6%. The particle identification effi¬ 
ciency is studied with the sample J/ip —>■ (pp —>■ 
K+K- ■ 77 . The efficiency difference between data 
and MC is 0.5%. To estimate the uncertainty asso¬ 
ciated with the kinematic constraint, a control sam¬ 
ple of J/ip (j)r] ^ Ar+ 7 f “77 is selected without a 
kinematic fit. The efficiency is the ratio of the sig¬ 
nal yields with and without the kinematic require¬ 
ment x^(4C') < 30. The difference between data 
and MC, 3.2%, is assigned as the systematic uncer¬ 
tainty. For the uncertainties from the fit, alternative 


fits are performed by varying the bin size and fit 
ranges. In addition, we also consider the effect from 
the parameterization of the function (Fhq) for the 
null hypothesis, the relative phase angle S and the 
decay width r(m). We repeat hts parameterizing 
Fho with a third-order polynomial and extending S 
in Eq. nto be <5 ± -b , 52 with two 

more parameters and 62 - Assuming the modes 
(f) —>■ K~^ K~ /KsKl have the same branching frac¬ 
tion 50%, we also perform a fit replacing r(m) with 
Tx+im) X 50% + rjio{m) x 50%, where rK+/Ko{m) 
is Eq. [3] using the mass of K'^/K^. The yield differ¬ 
ence with respect to the nominal fit is taken as the 
systematic uncertainty due to the parameterization. 
The mass resolution, = I.00± 0.02 MeV/c^, is 
determined from J/ip —> Varying cr^ within 

±0.02 MeV/c^ in the fit, the signal yield difference 
compared to the nominal fit is less than 1%. The 
QED contribution is neglected and the uncertainty 
for Ai^°o (3.097) is taken as 5.8 as stated above. It 


contributes a systematic uncertainty of 0.7% (16.4%) 
for Solution I (II), ignoring the possible interference 
between the QED process and J/ip resonance decay. 
The mass and width of the (j) meson have been hxed 
to their world averages Changing them with 

1(7 uncertainty, the signal yield difference is taken 
as the systematic uncertainty. The meson radius r 
is 3 GeV“^ in the nominal fit. We change it from 
1 GeV“^ to 5 GeV“^, and the largest signal yield 
difference is 2.3% (3.0%) for Solution I (II). In the 
ht, the 7 r° peaking background J/'tj) —>■ ^r]c{lS) —>■ 
'yK'^K~TY^ is neglected. These background events 
can be subtracted by a MC simulation normalized 
according to the relevant branching fractions 17 1 
and the efficiency. The signal yield difference is 1.1% 
(3.8%) for Solution I (II). We also consider the un¬ 
certainties from the number of J/'ip events and the 
branching fraction of ^ K^K~. The total sys¬ 
tematic uncertainty in Table Hill is the quadratic sum 
of the individual ones, assuming they are indepen¬ 
dent. 


In summary, based on 1.31 billion J/'i/ events col¬ 
lected with the BESIII detector, we perform an anal¬ 
ysis of the decay J/ijj —>■ —)> and find 

a structure around 1.02 GeV/c^ in the K'^K~ in¬ 
variant mass spectrum. It can be interpreted as 
interference of J/(/> —>■ (/ 7 r° with other processes 
decaying to the same final state. The fit yields 
two possible solutions and thus two branching frac¬ 
tions, [2.94 ± 0.16(stat.) ± 0.16(syst.)] x 10“® and 
[1.24±0.33(stat.)±0.30(syst.)] x 10"^. 

Ref. 0 provides a model-independent 
relation, = (r_Etan 0 y — 

l/V2)V(tan6»y/\^ ± rs)'^. Here B{VP) = 
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TABLE III. Summary of branching fraction systematic 
uncertainties (in %). 


Source 

Solution I Solution H 

MDC tracking 

2.0 

2.0 

Photon detection 

2.0 

2.0 

Particle identification 

0.5 

0.5 

7r°’s selection 

0.6 

0.6 

Kinematic fit 

3.2 

3.2 

Bin size 

1.0 

6.5 

Fit range 

1.0 

15.3 

Mass resolution 

0.1 

0.4 

Parameterization of Fhq 

0 

1.9 

Parameterization of S 

0.9 

1.6 

Parameterization of r(m) 

0.1 

0.0 

QED continuum 

0.7 

16.4 

The mass and width of (p 

0.8 

0.1 

The meson radius r 

2.3 

3.0 

7r° peaking background 

1.1 

3.8 

Number of J/ip 

0.8 

0.8 

Uncertainty of B{(p —>■ K'^K~] 

1 1.0 

1.0 

Total 

5.5 

24.4 



Te Qy (degree) 

FIG. 3. Dependence of the reduced branching fraction 
ratio (a) on the nonet symmetry break¬ 

ing strength te assuming lj — (j) ideal mixing and (b) 
on the mixing angle 6v assuming nonet symmetry. The 
yellow (green) box represents Solution I (II). The blue 
line represents the nonet symmetry value in (a) and the 
ideal mixing angle in (b). 


B{VP)/py is the reduced branching fraction of 
the decay J/tp VP, and pv is the momentum 
of the vector meson in the rest frame of J/ip', Oy 
is the (jj — (j) mixing angle; te is a dimensionless 


parameter accounting for nonet symmetry break¬ 
ing in the electromagnetic sector and te = 1 
corresponds to nonet symmetry. We have used 
wttO) = (4.5 ± 0.5) X [l^. li uj - p 

are mixed ideally, namely Oy = = arctan , 

the nonet symmetry breaking strength is 
6 e = te-I = (+21.0 ± 1.6)% or (-16.4 ± 1.0)% 
((+3.9 ± 0.8)% or (-3.7 ± 0.7)%) for Solution I 
(II), illustrated in Fig. [3] (a). On the other hand, 
we obtain ^y = \9y - = 4.97° ± 0.33° 

(1.03° ± 0.19°) for Solution I (II) assuming nonet 
symmetry, shown in Fig. |3] (b). However, py is 
found to be 3.84° from the quadratic mass formu¬ 
lae [13 and 3.34° ± 0.09° from a global fit to the 
radiative transitions of light mesons [l^]- The (py 
values do not agree with either solution. This is the 
first indication that nonet symmetry [3 is broken 
and the doubly OZI-suppression process contributes 
in J/ip electromagnetic decays. 
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